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Abstract 

We compute in linearized gravity all the contributions to the gravitational effective action due 
to a virtual Dirac fermion, related to the conformal anomaly. This requires, in perturbation theory, 
the identification of the gauge-gauge-graviton vertex off mass shell, involving the correlator of the 
energy-momentum tensor and two vector currents (TJJ), which is responsible for the generation 
of the gauge contributions to the conformal anomaly in gravity. We also present the anomalous 
effective action in the inverse mass of the fermion as in the Eulcr-Heisenberg case. 
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1 Introduction 



Investigations of conformal anomalies in gravity (see [T] for an historical overview and references) [2] 
and in gauge theories [31 SI [5] as well as in string theory, have been of remarkable significance along 
the years. In cosmology, for instance, [6] (see also [7] for an overview) the study of the gravitational 
trace anomaly has been performed in an attempt to solve the problem of the "graceful exit" (see for 
instance [SJ [OJ [TU| lll|). In other analysis it has been pointed out that the conformal anomaly may 
prevent the future singularity occurrence in various dark energy models [UJ [13] . 

In the past the analysis of the formal structure of the effective action for gravity in four dimen- 
sions, obtained by integration of the trace anomaly [141 H5| . has received a special attention, showing 
that the variational solution of the anomaly equation, which is non-local, can be made local by the 
introduction of extra scalar fields. The gauge contributions to these anomalies are identified at 1-loop 
level from a set of diagrams - involving fermion loops with two external gauge lines and one graviton 
line - and are characterized, as shown recently by Giannotti and Mottola in [16], by the presence 
of anomaly poles. Anomaly poles are familiar from the study of the chiral anomaly in gauge theo- 
ries and describe the non-local structure of the effective action. In the case of global anomalies, as in 
QCD chiral dynamics, they signal the presence of a non-perturbative phase of the fundamental theory, 
with composite degrees of freedom (pions) which offer an equivalent description of the fundamental 
lagrangian, matching the anomaly, in agreement with 't Hooft's principle. Previous studies of the 
role of the conformal anomaly in cosmology concerning the production of massless gauge particles and 
the identification of the infrared anomaly pole are those of Dolgov [171 [T8] , while a discussion of the 
infrared pole from a dispersive derivation is contained in [19] . 

In a related work [20] we have shown that anomaly poles are typical of the perturbative description 
of the chiral anomaly not just in some special kinematical conditions, for instance in the collinear 
region, where the coupling of the anomalous gauge current to two (on-shell) vector currents (for the 
AVV diagram) involves a pseudoscalar intermediate state (with a collinear and massless fermion- 
antifermion pair) but under any kinematical conditions. They are the most direct - and probably also 
the most significant - manifestation of the anomaly in the perturbative diagrammatic expansion of 
the effective action. On a more speculative side, the interpretation of the pole in terms of composite 
degrees of freedom could probably have direct physical implications, including the condensation of the 
composite fields, very much like Bose Einstein (BE) condensation of the pion field, under the action 
of gravity. Interestingly, in a recent paper, Sikivie and Yang have pointed out that Peccei-Quinn 
axions (PQ) may form BE condensates [21J. With these motivations in mind, in this work, which 
parallels a previous investigation of the chiral gauge anomaly [20], we study the perturbative structure 
of the off-shell effective action showing the appearance of similar singularities under general kinematic 
conditions. Our investigation is a first step towards the computation of the exact effective action 
describing the coupling of the Standard Model to gravity via the conformal anomaly, that we hope to 
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discuss in the future. 

In our study we follow closely the work of |16j . There the authors have presented a complete off- 
shell classification of the invariant amplitudes of the relevant correlator responsible for the conformal 
anomaly, which involves the energy momentum tensor (T) and two vector currents (J), TJJ, and 
have thoroughly investigated it in the QED case, drawing on the analogy with the case of the chiral 
anomaly. The analysis of [TB] is based on the use of dispersion relations, which are sufficient to identify 
the anomaly poles of the amplitude from the spectral density of this correlator, but not to characterize 
completely the off-shell effective action of the theory and the remaining non-conformal contributions, 
which will be discussed in this paper. The poles that we extract from the complete effective action 
include both the usual poles derived from the spectral analysis of the diagrams, which are coupled in 
the infrared (IR) and other extra poles which account for the anomaly but are decoupled in the same 
limit. These extra poles appear under general kinematic configurations and are typical of the off-shell 
as well as of the on-shell effective action, both for massive and massless fermions. 

We also show, in agreement with those analysis, that the pole terms which contribute to the 
conformal anomaly are indeed only obtained in the on-shell limit of the external gauge lines, and 
identify all the mass corrections to the correlator in the general case. This analysis is obtained 
by working out all the relevant kinematical limits of the perturbative corrections. We present the 
complete anomalous off-shell effective action describing the interaction of gravity with the photons, 
written in a form in which we separate the non-local contribution due to the anomaly pole from the 
rest of the action (those which are conformally invariant in the massless fermion limit). Away from 
the conformal limit of the theory we present a 1/m expansion of the effective action as in the Euler- 
Heisenberg approach. This expansion, naturally, does not convey the presence of non-localities in the 
effective action due to the appearance of massless poles. 

The computation of similar diagrams, for the on-shell photon case, appears in older contributions 
by Berends and Gastmans |22j using dimensional regularization, in their study of the gravitational 
scattering of photons and by Milton using Schwinger's methods [23] , The presence of an anomaly 
pole in the amplitude has not been investigated nor noticed in these previous analysis, since they 
do not appear explicitly in their results, nor the 1/m expansion of the three form factors of the on- 
shell vertex, contained in [22], allows their identification in the S-matrix elements of the theory. Two 
related analysis by Drummond and Hathrell in their investigation of the gravitational contribution to 
the self-energy of the photon [24J and the renormalization of the trace anomaly |25j included the same 
on-shell vertex. Later, this same vertex has provided the ground for several elaborations concerning a 
possible superluminal behaviour of the photon in the presence of an external gravitational field [26j . 
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2 The conformal anomaly and gravity 



In this section we briefly summarize some basic and well known aspects of the trace anomaly in 
quantum gravity and, in particular, the identification of the non-local action whose variation generates 
a given trace anomaly. 

We recall that the gravitational trace anomaly in 4 spacetime dimensions generated by quantum 
effects in a classical gravitational and electromagnetic background is given by the expression 

1 
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2bC 2 + 2b' E - -OR +2cF 



(1) 



where b, b' and c are parameters that for a single fermion in the theory result b = 1/320 ir 2 , b' = 
— 11/5760 7r 2 , and c = — e 2 /247r 2 ; furthermore C 2 denotes the Weyl tensor squared and E is the Euler 
density given by 

C 2 = C^ vp C x ^ = R^ up R x ^-2R^ + ^- (2) 

E = *R\ t j,up *R XtJ,vp = R\ fJiUp R XfJ '' /p - ^RfiuR' 1 " + R 2 - (3) 

The effective action is identified by solving the following variational equation by inspection 

2 5T 

Its solution is well known and is given by the non-local expression 



Sanomidt A] — (5) 

i J (fixy/=g~ j (fix'sf^g 1 (e - ^DRj G 4 (x,x') 2b C 2 + b' (e - -DRJ + 2cF flu F flu . 

Notice that we are omitting ^fgR 2 terms which are not necessary at one loop level. The notation 
G^XjX 1 ) denotes the Green's function of the differential operator defined by 



A 4 = V M (v»V u + 2R tMU - = D 2 + 2R^V fl V u + -(V^V,, - -BD 



(6) 



and requires some boundary conditions to be specified. This operator is conformally covariant, in fact 
under a rescaling of the metric one can show that 



9nu = e a g„ v -»■ A 4 = e~ Z(7 A 4 . (7) 

Notice that the general solution of © involves, in principle, also a conformally invariant part that is 
not identified by this method. As in ref. |16j . we concentrate on the contribution proportional to F 2 
and perform an expansion of this term for a weak gravitational field and drop from this action all the 
terms which are at least quadratic in the deviation of the metric from flat space 

9pu =V(iv + K h pu K 2 = 16 7T G, (8) 
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with G the gravitational constant. The non-local action reduces to 

Sanom[9,A] =~J rf 4 ^ J d 4 x'y^R^ [F af) F a %, , (9) 

valid for a weak gravitational field. In this case 

R ( x ] = ^ dl hT -ah, h = h^. (io) 

The presence of the Green's function of the □ operator in Eq. is the clear indication that the 
solution of the anomaly equation is characterized by an anomaly pole. In the next sections we are 
going to perform a direct diagrammatic computation of this action and reobtain from it the pole 
contribution identified in the dispersive analysis of [16] and the conformal invariant extra terms which 
are not present in ([9]). We start with an analysis of the correlator following an approach which is 
close to that followed in ref. [16]. The crucial point of the derivation presented in that work is the 
imposition of the Ward identity for the TJJ correlator (see Eq. ([42 p below) which allows to eliminate 
all the Schwinger (gradients) terms which otherwise plague any derivation based on the canonical 
formalism and are generated by the equal-time commutator of the energy momentum tensor with the 
vector currents. In reality, this approach can be bypassed by just imposing at a diagrammatic level the 
validity of an operatorial relation for the trace anomaly, evaluated at a nonzero momentum transfer, 
together with the conservation of the vector currents on the other two vector vertices of the correlator. 

3 The construction of the full amplitude r^ val3 (p, q) 

We consider the standard QED lagrangian 



C = ~F tlu F^+i^(d t t-ieA^-mW, (11) 



with the energy momentum tensor split into the free fermionic part Tf, the interacting fermion-photon 
part Tf p and the photon contribution T p h which are given by 

Tf = -iijr/^d^ + g^ u {iipj x d\^ - mi/nl>), ( 12 ) 
Tf = - eJ^A u ^ + egT J X A X , (13) 



and 



where the current is defined as 



J*(x) = ^{x)^^{x) . (15) 
In the coupling to gravity of the total energy momentum tensor 
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we keep terms linear in the gravitational field, of the form h^T^, and we have introduced some 

standard notation for the symmetrization of the tensor indices and left-right derivatives H^ u ) = 

•o- -» <— 

(H^ v + H Ufl )/2 and d^ = {d^— d/J/2. It is also convenient to introduce a partial energy momentum 
tensor T p , corresponding to the sum of the Dirac and interaction terms 

Tg v = Tf + T% (17) 

which satisfies the inhomogeneous equation 

ajr? = -dji%. (is) 

Using the equations of motion for the e.m. field d u F^ u = J 1 *, the inhomogeneous equation becomes 

d v Tf = F» X J X . (19) 

There are two ways to identify the contributions of T^ u and Tp V in the perturbative expansion of the 
effective action. In the formalism of the background fields, the T p JJ correlator can be extracted from 
the defining functional integral 

{Tf(z)) A = J D^jDi) T£ u {z) e i!d 4 xC+fJ-A{x)d±x 

= (Tg v e i f d4xJ - A (^ (20) 

expanded through second order in the external field A. The relevant terms in this expansion are 
explicitly given by 

PTWU = ^{z){j -a)(j. a)) + (t;;(j • a)> + ... , (21) 

with (J ■ A) = J d A xJ ■ A(x). The corresponding diagrams are extracted via two functional derivatives 
respect to the background field A^ and are given by 



- 8A a { X )5A p {y) 



A=0 



V^afS = ( - e)2 (Ty{z)J a (x)jP{y)) A=0 (23) 



W raf> = ^{z){J-A)) 



8A a (x)8Ap(y) 



A=0 



= 5\x - z)g a ^Il u) P(z, y) + 5\y - z)g^U u)a {z, x) - g^[6 4 {x - z) - 5 4 (y - z)]U a ^{x, y) 

(24) 

These two contributions are of 0(e 2 ). Alternatively, one can directly compute the matrix element 

MT = {0\Tf(z) J d 4 wd 4 w'J ■ A(w)J • A(u/)|77>, (25) 
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which generates the diagrams (b) and (c) shown in FigfTJ respectively called the "triangle" and the 
"t-p-bubble" ("t-" stays for tensor), together with the two ones obtained for the exchange of p with q 
and a with j3. 

The conformal anomaly appears in the perturbative expansion of T p and involves these four di- 
agrams. The electromagnetic contribution is responsible for other two diagrams whose invariant 
amplitudes are well-defined and will be used to fix the ill-defined amplitudes present in the tensor 
expansion of Tp U by using a Ward identity. 

The lowest order contribution is obtained, in the background field formalism, from Maxwell's e.m. 
tensor, and is given by 



SA a (x)5Ap(y) 



(26) 

A=0 



Equivalently, it can be obtained from the matrix element 

<o|2yiry> (27) 

which allows to identify the vertex in momentum space. Using (|26p we easily obtain 



S[j,vaf3(z,x,y) = 2g a pd( yfJL 5 xz d u )5 yz — 2g^ fl d u ^S xz d a 5 yz — 2g a ^ v d ^8 yz d^8 xz 

+9ati9/3udxS yz d 5 XZ + gavQp^dx^yzd 5 XZ + g^pdpdxzdadyz — d p 5 yz d p 

(28) 

where <9^<5 X2 = d/dx^5{x — z) and so on. In momentum space this lowest order vertex is given by 
gi^afi = (^/y + g°0 +p . q (g a » ' g^ + g a »g^) - g^ (p ■ q g a/3 - q a p P ) 

- (/V + /V) g a - {g^q 11 + g afl q u ) (29) 

The corresponding vertices which appear respectively in the triangle diagram and on the t-bubble 
at 0(e 2 ) are given by 

V'^(h,k 2 ) = ±[^(k l + k 2 y+Y(ki + k 2 r]-^[ 1 x (k 1 + k 2 ) x -2m}, (30) 
w'^ a = -i( 7 V a + 7V a ) + <T7", (31) 

where k\{k 2 ) is outcoming (incoming). Using the two vertices V ,fMU (ki, k 2 ) and W'^ ua we obtain for 
the diagrams (b) and (c) of FigQ] 



V^(p,q) = -(-iefi 3 



dH tr {V'^(l +p, l-q)(f,-i + m)^ (jt + m) j a (J + $ + m)} 



(2vr) 4 [I 2 - m 2 } [(I - q) 2 - m 2 ] [(I + p) 2 - m 2 } 

(32) 
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Figure 1 : The complete one-loop vertex (a) given by the sum of the 1PI contributions called V^ va " (p, q) (b) 
and W^ al3 { Pl q) (c). 



and 

T ™/3, v ,.2 V 2 f dH tr{W'^ a (} + m) 7 ^-4 + m)} 

W M = ~ {le)l 7(2^ p2- m 2 ][a _ g) 2_ m 2] > ^ 

so that the one-loop amplitude in Fig. [T] results 

(p, q) = (p,q)+ V^ a (q,p) + (p,q)+ W^ a (q,p). (34) 

The bare Ward identity which allows to define the divergent amplitudes that contribute to the anomaly 
in r in terms of the remaining finite ones is obtained by re-expressing the classical equation 

dvTft = -F^J U (35) 

as an equation of generating functionals in the background electromagnetic field 

du{T^) A = -F^(J V ) A , (36) 

which can be expanded perturbatively as 

d v (Tg) A = -F^{J V J d 4 w{ie)J ■ A(w))+ ... . (37) 

Notice that we have omitted the first term in the Dyson's series of {J V ) A , shown on the r.h.s of ([3 
since (J v ) = 0. The bare Ward identity then takes the form 



5" (F^(z)(J x (z)) A ) 



6A a (x)6Ap(y) 



(38) 

A=0 



which takes contribution only from the first term on the r.h.s of Eq. (|37|) . This relation can be written 
in momentum space. For this we use the definition of the vacuum polarization 

U^(x,y) = -ie 2 (J a (x)My)), (39) 
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p^pVq a pfi 


qU jpQL q(3 


qfipU q(y.pj3 


qVpVq a qP 


gOLfl gfiv 




pVqVpCpP 


q^p v p a qP 


qfi qV p/3 


q^q u p a qP 


gOtV gfiH 




qfi pis p/3 






q^q 1 'q a p 13 




p»p v g a P 


pPp v g a » 


pPpVg av 


p a p v g^ 


pVpagPv 


P» p 0g^ 






pPqftgOv 


p a q u g^ 


pHq* g P» 


p a q f3 g^ v 


q^p v g a P 


q P p Vg<XH 


q Pp(t g ov 


qa p VgPv 


qVp* g P» 


q a pPg^ v 


q»q v g af3 


qpqVgay, 


qPq^g au 


q a q u g^ 


q fl q a gP u 


q a qPg^ 







Table 1 : The 43 tensor monomials built up from the metric tensor and the two independent momenta p and q 
into which a general fourth rank tensor can be expanded. 



or 



n a/3 ( P ) 



-i 2 {-ief 



dH tr{ 7 a (/ + m)7 /3 (/+^ + m)} 
(2vr) 4 [I 2 - m 2 ] [(I + p) 2 



(p 2 g ap -p a p f} )Tl{p 2 ,m 2 ^ 



with 



U(p 2 ,m 2 ) 



QA (m 2 )+p 2 -Qm 2 - 3B Q (p 2 ,m 2 ) (2m 2 + p 2 ) 



(40) 



(41) 



36 7T 2 p 2 

which obviously satisfies the Ward identity p a H a ^(p) = 0. The expressions of the .4.0 an d Bq contri- 
butions are given in Appendix [XI 

Using these definitions, the unrenormalized Ward identity which allows to completely characterize 
the form of the correlator in momentum space becomes 



k„T^(p,q) 



fp a pP - q^g^p 2 + g^q a p 2 - g^p a p • q) n(p 2 ) 
+ (V - pVV + g m p p q 2 - g m <fp ■ g) i% 2 ) . 



(42) 



3.1 Tensor expansion and invariant amplitudes of T 



The full one-loop amplitude T can be expanded on the basis provided by the 43 monomial tensors 
listed in Tabffl 

43 

r^(p,q) = £ A(k 2 ,p 2 ,q 2 )ir^P,q)- (43) 



i=l 



Since the amplitude r fJ-ua ^(p,q) has total mass dimension equal to 2 it is obvious that not all the 
coefficients Ai are convergent. They can be divided into 3 groups: 
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a) A\ < Ai < Aiq - multiplied by a product of four momenta, they have mass dimension —2 and 
therefore are UV finite; 

b) An < Ai < A\g - these have mass dimension 2 since the four Lorentz indices of the amplitude 
are carried by two metric tensors 

c) A20 < Ai < A^ 3 - they appear next to a metric tensor and two momenta, have mass dimension 
and are divergent. 

The way in which the 43 invariant amplitudes will be managed in order to reduce them to the 13 
named Fi(k 2 ,p 2 ,q 2 ) is the subject of this section. The reduction is accomplished in 4 different steps 
and has as a guiding principle the elimination of the divergent amplitudes Ai in terms of the convergent 
ones after imposing some conditions on the whole amplitude. We require 

a) the symmetry in the two indices \i and v of the symmetric energy- momentum tensor T^ v \ 

b) the conservation of the two vector currents on p a and q@ ; 

c) the Ward identity on the vertex with the incoming momentum k defined above in Eq. (|3.ip . 
Condition a) becomes 

T ^ap (p :q ) = Y v ^ (p, q) , (44) 

giving a linear system of 43 equations; 15 of them being identically satisfied because the tensorial 
structures are already symmetric in the exchange of (i and u, while the remaining 14 conditions are 

^4i8 = A19, 
A 27 = A 31 , 

(45) 

where all Ai are thought as functions of the invariants k 2 ,p 2 ,q 2 . After substituting ([43]) into T fiua ^(p, q) 
the 43 invariant tensors of the decomposition are multiplied by only 29 invariant amplitudes. Condition 
b), which is vector current conservation on the two vertices with indices a and (3, allows to re-express 
some divergent Ai in terms of other finite ones 

p a (p, q) = qp (p, q) = 0. (46) 

This constraint generates two sets of 14 independent tensor structures each, so that in order to fulfill 
each coefficient is separately set to vanish. The first Ward identity leads to a linear system 



.45 


= Aq, 


As 


= A 9 , 


A 10 


= A u , 


A 13 


= A u 


A n 


= A 22 , 


A 2 4 


= A 28 , 


A 25 


= A 29 , 


A 26 


= A 30 


A 32 


= A 36 , 


-434 


= A 37 , 


A 33 


= A 38 , 


A 35 


= A 39 
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composed of 10 equations 



p a T^(p,q)=0 



(47) 



A 19 + A 36 p ■ p + A 37 p ■ q = 0, 
A 38 p-p + A 39 p- q = 0, 
A 17 + A 40 p ■ p + A i2 P ■ q = 0, 
A 44 p-p + A 43 p- q = 0, 
^20 + 2A 28 + Ai p ■ p + A 4 p ■ q = 0, 
2^3o + A 3 p-p + ^ 7 p-g = 0, 
A 2 2 + -429 + ^6 P ■ P + -Ail P ■ q = 0, 
A 31 + A 9 p-p + A 14 p ■ q = 0, 
^23 + -4i2 P ■ P + -4i6 P ■ 9 = 0, 

we choose to solve it for the set {A15, An, Ai 9 , A23, A28, A29, A 3 q,A 3 i,A 3 $, A 43 } in which only the first 
one is convergent and the others are UV divergent. The set would not include all the divergent Ai 
since in the last equations appear two convergent coefficients, A15 and A2. 
Following our choice the result is 



A 15 = -A 2 P -«, 
p ■ p 

-A 36 p- p - A 37 p - q 

1 - 

2 . 



A 19 
A 2S 
A 30 
A 39 



120 



A\p • p — A 4 p • q 



1 

2 L 



-A 



A 3 p-p + A 7 p- q 
p-p 



38" 



p-q 



An = ~A m p ■ p - A 42 p ■ q, 
A23 = -A12 p-p- A w p-q, 
A 29 = -A22 - A 6 p-p- Anp-q, 

A 31 

A 43 = -A 41 



-A 9 p ■ p - A 14 p ■ q, 
p-p 



p-q 



(48) 
(49) 
(50) 

(51) 
(52) 



In an analogous way we go on with the second Ward identity (WI) after replacing the solution of 
the previous system in the original amplitude. The new one is indicated by (p, q), where the 

subscript b is there to indicate that we have applied condition b) on T. The constraint gives 

A 40 p ■ q + A 41 q ■ q =0, 
Aip ■ q + A 3 q ■ q =0, 
A 20 + A 4 p-q + A 7 q-q =0, 
A 36 + A 6 p ■ q + A 9 q ■ q =0, 
^22 + A 37 + A u p ■ q + A 14 q ■ q 



q^r P ip,Q) 



(53) 



0, 



2A 38 + A 12 p-q-A2^=0. 



We solve these equations determining the amplitudes in the set {A±, A20, A22, A 3 q, A 38 , A 4 o} in terms 
of the remaining ones, obtaining 

Ayip-pp-q -A 2 p-q q-q A 4 \q-q 



A 3% = -'- 



2p ■ p 



A 



40 



p-q 



(54) 
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A-v 
A 22 



A 3 q-q 

p-q 

-A 37 - A u p-q -Auq-q. 



A 20 = -A 4 p ■ q - A 7 > 
A 3e = -A 6 p-q - A 9 c 



(55) 
(56) 



The manipulations above have allowed a reduction of the number of invariant amplitudes from the 
initial 43 to 13 using the {^t,^} symmetry (14 equations), the first WI on p a (10 equations) and the 
second WI on qp (6 equations). The surviving invariant amplitudes in which the amplitude Tc Ua ^(p, q) 
can be expanded using the form factors are {A2, A 3 , A 4 , Aq, A 7 , Aq, Ah, A12, A14, A±q, A 37 , A41, A^}- 
This set still contains 3 divergent amplitudes, (A 37 , A41, A43). The amplitude Fc ua ^ (p, q) is indeed 
ill-defined until we impose on it condition c), that is Eq. (142 p . This condition gives 



Eq.g2]) 



-A, 



1 + 



p-p 

2p-q 



+ A 6 + lA 7 -A 9 - If = 



A 37 + A 42 +A 4 \p-p+p-q]+A n p-q+±A 7 q-q + 

+ A n q-q+kA 3 Z%fL = 0, 
\A2^-A 4l ^±fL-\A 3 p.p + A 7 {p.p+\p.q)+A % p.q 



p-q 



+A l2 {\p ■ q + q-q)+ A U ( P -q + 2q-q) + 2A 37 - U(p 2 ) - U(q 2 ) = 



From this condition we obtain 
A 37 -- 



A 2 p-qq-q 1 1 . \ t 1 a fP'P + VQ 

— V-A 3 p-p - ~A 7 {2p-p + p-q) + -A 41 [ 

4p-p 4 4 2 \ p ■ q 

■~A 6 p-q- -A l2 (p-q + 2q-q)- \a 14 (p-q + 2q-q) + \ [U(p 2 ) + U(q 2 ) 



A?, 

I41 = —j-P-P- {A 3 -A 6 -A 7 + A 9 )p-q 



A42 = -J- P ■ P 



p • p 



+ 1 



+ 5 A 7 (p ■ p + p ■ q - q ■ q) - A 4 (p-p + p 



p-q v 

-(Aj - A 9 )p ■ p + (A u - A n )(q -q + p 



(57) 
(58) 

(59) 



After these steps we end up with an expression for V written in terms of only 10 invariant amplitudes, 
that are X = {j4 2 , j4 3 , j4 4 , A> A> ^9) ^lii ^leli significantly reduced respect to the original 

43. Further reductions are possible (down to 8 independent invariant amplitudes), however, since 
these reductions just add to the complexity of the related tensor structures, it is convenient to select 
an appropriate set of reducible (but finite) components characterized by a simpler tensor structure 
and present the result in that form. The 13 amplitudes introduced in the final decomposition are, in 
this respect, a good choice since the corresponding tensor structures are rather simple. These tensors 
are combinations of the 43 monomials listed in TabfTJ 

The set X is very useful for the actual computation of the tensor integrals and for the study of 
their reduction to scalar form. To compare with the previous study of Giannotti and Mottola |16| we 
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have mapped the computation of the components of the set X into their structures Fi (i = 1, 2, .., 13). 
Also in this case, the truly independent amplitudes are 8. One can extract, out of the 13 reducible 
amplitudes, a consistent subset of 8 invariant amplitudes. The remaining amplitudes in the 13 tensor 
structures are, in principle, obtainable from this subset. 

3.2 Reorganization of the amplitude 

Before obtaining the mapping between the amplitudes in X and the structures Fi, we briefly describe 
the tensor decomposition introduced in [16J which defines these 13 structures. We define the rank-2 
tensors 

u a ^p,q) = {p-q)g^-q cx p^ (60) 
w a/3 (p, q) = p 2 q 2 g a/3 + {p-q) p a q P - q 2 p a / - p 2 q a q P , (61) 

which are Bose symmetric, 

u^(p,q)= U ^(q,p), (62) 

w^(p,q)=w^(q,p), (63) 

and conserve vector current, 

PaU a P(p,q) =qpu a P(p,q) =0, (64) 

Pa w a ? (p, q) = qp w a/3 (p, q) = . (65) 

These two tensors are used to build the set of 13 tensors catalogued in Table [5J They are linearly 
independent for generic k 2 ,p 2 ,q 2 different from zero. Five of the 13 tensors are Bose symmetric, 
namely, 

t^ipA)=t^ a {q,p), i = 1,2,7,8,13, (66) 
while the remaining eight tensors form four pairs which are overall related by Bose symmetry 

C /, (P.9)=C /to (9.P). (67) 
tr^P,q)=t^ a (q,p), (68) 
t^( P ,q)=t^ a (q,p), (69) 

tr*(p.?)=*&*w)- (™) 

The amplitude in (j34"|) can be expanded in this basis composed as 

13 

T^( P ,q) = J>( S ; Sl)S2 ,m 2 ) t^(p,q), (71) 

8=1 
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i 




1 


(fcV - A;^) u Q/3 (p.g) 


2 


{k 2 g^ v - Wk u ) w al3 {p.q) 


3 


(p 2 g ^ - Ap»p u ) u af3 {p.q) 


4 


[p 2 g^ -4p»p u ) w al3 {p.q) 


5 


(q 2 g^ -4q»q u ) u a0 (p.q) 


6 


(q 2 gf MU - AqVq v ) w a/3 {p.q) 


7 


[p -qsT- 2[<fp v + p»q v )\ u a/3 (p.q) 


8 


[p -qg^ - 2(g*y + jfq")] w a P{p.q) 


9 


(p.qp a - p 2 q a ) \pP (q^p v + p^q u ) -p-q (gP v p» + g^p u )} 


10 


{p-qqP - q 2 p P ) [q a {q^p u + p^q v ) -p-q {g au q^ + g afM q u )] 


11 


(p ■ qp a - p 2 q a ) [2 qPq^q" - q 2 {g^q^ + g^q v )\ 


12 


(p-qq 13 - gV) [2p a p fl p u - p 2 (g av p^ + g a ^p v )] 


13 


(p)* q v + p» q »)g<*P +p . q {g av g^ + ~ ^"u^ 
-(gPvj? + g^p u )q a - {g av q^ + g a »q v )pP 



Table 2: Basis of 13 fourth rank tensors satisfying the vector current conservation on the external 
lines with momenta p and q. 



where the invariant amplitudes Fi are functions of the kinematical invariants s = k 2 = (p + q) 2 , 
si = p 2 , S2 = q 2 and of the internal mass m. In [16] the authors use the Feynman parameterization 
and momentum shifts in order to identify the expressions of these amplitudes in terms of parametric 
integrals, which was the approach followed also by Rosenberg in his original identification of the 6 
invariant amplitudes of the AVV anomaly diagram. If we choose to reorganize all the monomials into 
the simpler set of 13 tensor groups shown in TabJJl then we need to map the A{ in \ an d the F^s. 
The mapping is given by 

1 



3 k 2 



At(4p- q + 3p- p) + 2A u (p-q + 2q-q) + 2A 6 p-p 

p-p q-q 



+ 2A 7 q-q - 2A 14 q ■ q - A 16 q ■ q + 2A 3 



F 2 
F 3 

Fa 



1 



3 k 2 
1 

12 A; 2 



-2A, 



p ■ p 
p-q 



+ 2 + 4A 6 + A 7 



p-q 

2A 9 - A12 



[A 4 (2p-q + 3q-q)-2A ll {p-q + 2 q ■ q ) - 2A 6 p ■ p 



2A 7 q ■ q + 2A 14 q-q + A 16 q-q-2A 3 
p ■ p 



A 7 1 

+ 



4p-p 12 k 2 



2A, 



p-q 



+ 2 



4A f 



p ■ p q ■ q 
p-q 

-A 7 + 2A 9 + Au 



(72) 
(73) 

(74) 
(75) 
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Aw 1 
4 12 k 2 



-2A G p-p -2A 3 Q ' qP ' P +A 4 (-3p-p -4p 
p-q 



-2 A n (p ■ q +2q ■ q) - 2 A 7 q ■ q + 2 A 14 q ■ q + A w q ■ q] , 



F 6 
F 7 



A l2 1 



-4A 6 -A 7 + 2A 9 + A 12 + 2A 3 ( ^+ 2 

p-q 



Aq-q 12 k 2 

A n 1 ( . A 6 A u \ 1 

— + J^{ A9q ' qp ' p + T p ' pp ' q+ — q ' qp ' q ) + 6k 2 l 



(76) 
(77) 



A^-ip-q -3p-p] 



2A u (p-q + 2q-q)-2A 6 p-p - 2 A 7 q ■ q + 2A u q-q + A 16 q-q - 2 A 



p ■ p 



p-q 



F 8 
F 9 
F\o 
F u 

F\2 

F13 



1 



6 k 2 
Ae_ 
p-q 



2 A-, 



Aa 



P-P 

p-q 

q ■ q 

p-q 2 



+ 2 



Ag 
p-q 



{p-p + q-q)- 4:A 6 - A 7 - 4A 9 + A 



12 



. P-P . A u 

Ag - H 

p ■ q A p-q 
A12 A 2 



2 q ■ q 2p - p 
A 3 A 7 



2p ■ q 2p ■ p 
1 



(78) 
(79) 

(80) 

(81) 

(82) 

(83) 



A 2 p ■ qq ■ q 



Aq (p-p + p-q -q-q) + -A 7 (p-p+p-q-q-q) + 
2 4 ip ■ p 

+ -A 12 {p-q + 2q-q) + -^- 
4 4p • q 

+ -A g q-q + p-p h 1 

2 [ \p-q 



A 



ii 



p-q 



+ q-q 



[p-p 2 + (p- q +q ■ q)p- p +2p ■ qq ■ q) 



-[U(p)+U(q)]. 



(84) 



We have shown how to obtain the 13 Ff s, starting from our derivation of the one-loop full amplitude 
T IJ,ua ^(p, q) leading to the ten invariant amplitudes of the set X. Since we know the analytical expres- 
sion of the Ai involved, we can go one step further and give all the F£ s in their analytical form in 
the most general kinematical configuration. 



4 Trace condition in the non-conformal case 

Similarly to the chiral case, we can fix the correlator by requiring the validity of a trace condition 
on the amplitude, besides the two Ward identities on the conserved vector currents and the Bose 
symmetry in their indices. This approach is alternative to the imposition of the Ward identity (|42p 
but nevertheless equivalent to it. At a diagrammatic level we obtain 

g^T^ip, q) = A^(p, q) - ^ u^(p, q). (85) 
We comment below on this equation, in relation to the scales present in the perturbative expansion 
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of the correlator, which are, besides the fermion mass m, the energy at which we probe the correlator 
(s) and the subtraction point after renormalization (/j or M). We have also defined 



A a/3 (p, q) 



-m (ie) 2 J d A xd 4 ye ip - x+iq - y {^J a {x)J p (y)) 
dH 



-771 e 



tr 



a 1 1 

7 » T 



(27r) 4 + j> — m fi — m jt — </[ — m 
A direct computation gives 

A^(p, q) = G x {s, Sl ,s 2 ,m 2 ) u a P(p, q) + G 2 (s, Sl ,s 2 ,m 2 ) w a ?(p, q), 



+ exch. 



(86) 



(87) 



where 



3sFi(s,si,S2,m 2 ) 
3s F2(s,si,S2,m 2 ) 



Gi(s,si,s 2 ,m 2 
G 2 (s,si,s 2 ,m 2 



6vr 2 



and 



n ( 2n e l m . e i p2(s,S2,m 2 )s 2 m 2 f 2 .2,2,, i 

Gi(s,si,s 2 ,m J = — ^ 1 ^5 s + 4sis - 2s 2 s - 5si + s 2 + 4sis 2 

e 2 Pi(s,si,m 2 )sim 2 f 2 2 
s-n -(s-si) +5s 2 -4(s + si 



e 2 C (s,si,s 2 ,m 2 ) 



7T 2 <7 



s 2 J 

[(a - si) 3 - si + (3s + si) s 2 . + (-3s 2 - lOsis + s\) s 2 ] 
[(s-si) 2 -2s^ + (s + si)s 2 ] 



2m 7 1 



_2tt 2 ct 2 

2e 2 m 2 2e 2 P 2 (s, s 2 , m 2 )m 2 



ir 2 a 



(90) 



ir 2 a 2 



G 2 (s,si,s 2 ,m 2 ) = 

2e 2 Pi(s,si,m 2 )m 2 f 2 , 02,2, 
5—5 s + (si - 2s 2 ) s - 2sf + 4 + sis 2 

-e 2 Co(s,si,s 2 ,m 



H — [s - (Si + S 2 J S - [S 1 - 6s 2 Si + S 2 J S 



+ (si - s 2 ) 2 (si + s 2 ) 



(91) 



where 7 = s — si — S2, cr = s 2 — 2(si + s 2 ) s + ( s i ~~ s 2) 2 and the scalar integrals £>i(s, si, m 2 ), 
2?2(s,si,m 2 ), Co(s, si, S2, m 2 ) for generic virtualities and masses are defined in Appendix lAl 
We have checked that the final expressions of the form factors in the most general case, obtained either 
by imposing this condition on the energy momentum tensor or the Ward identity in the form given 
by Eq. (I36p exactly coincide. In Appendix |B] we discuss this relation in the simpler case of a massless 
fermion in the loop. 
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5 The off-shell massive (TJJ) correlator 



To obtain the explicit expression of the parametric integrals which describe the form factors, we follow 
an approach similar to that of |20] , for the case of the chiral gauge anomaly. These have been obtained 
by re-computing the anomaly diagrams by dimensional reduction together with the tensor-to-scalar 
decomposition of the Feynman amplitudes. For instance, in |20| we have given the explicit expressions 
of the parametric integrals of Rosenberg using this trick. The correctness of the result can be checked 
numerically by comparing the parametric forms to the explicit computation. In this case the procedure 
is identical, though the computations are very involved. By comparing the two approaches we extract, 
indirectly, an explicit expression of the parametric forms of these integrals, introduced in [16j . We have 
checked that indeed there is perfect numerical agreement between our computation and the parametric 
result, as discussed in Appendix O 

We introduce in this section the main results of our computation which will be used in the next 
sections for further analysis. The complete expressions of the form factors F{ (i = 1, . . . , 13) in the 
massive and then in the massless case are contained in Appendix iDl and [El respectively, whereas the 
master integrals are collected in Appendix [XJ In both cases the virtualities of the external lines are 
generic and denoted by s±, S2- After presenting the complete expressions, we discuss several kinematical 
limits of the result, in particular the on-shell limit for the two vector lines (si — >• 0, S2 — >■ 0) in order 
to better understand the structure of the whole correlator. The appearance of generalized anomaly 
poles in the correlator and their IR decoupling under the most general conditions will be discussed 
thoroughly. 

Notice that F13 contains two vacuum polarization diagrams with the two photon momenta which 
are divergent and we are bond to define a suitable renormalization of the 2-point function which will 
affect the running of the coupling. In the next section we will address the explicit relation between 
renormalization schemes and running of the coupling in the context of the renormalization of the 
correlator. 

5.1 Anomaly poles and their UV/IR significance 

There are close similarities between the effective action in the case of the chiral gauge anomaly and 
the conformal case, due to the presence of massless poles. In [20J we have analyzed the fact that 
in the chiral case the anomaly is entirely generated by the longitudinal component wl, which is 
indeed isolated for any configuration of the photon momenta. This is somehow unexpected since the 
dispersive analysis shows that the pole in wl is coupled only under a specific kinematic condition, 
and is usually interpreted as an infrared effect. Nevertheless there is a complete equivalence between 
the representation of the anomaly diagram in the Rosenberg representation - where the pole is not 
extracted as an independent component - and the L/T representation in which the pole is isolated 
under any kinematical configuration (and even in the massive case). This is apparent from the broken 
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anomalous Ward identities satisfied by the AVV diagram where the mass corrections and the anomaly 
term can be separately identified [20J. 

To illustrate the emergence of a similar behaviour in the case of the conformal anomaly, it is 
sufficient to notice in the expression of F± given in Eq. (|193p the presence of the isolated contribution 
(Fipoie = — e 2 /(187r 2 s)) which survives in the massless limit but is present also in the massive case. 
This component, indeed, is responsible for the trace anomaly also in the massive case, even though 
there appear extra corrections with mass-dependent terms. Obviously also in this case, which is generic 
from the kinematical point of view, one can clearly show that the pole does not couple in the infrared 
if we compute the residue of the entire amplitude. The anomaly pole, in fact, appears in the spectral 
function only in a special kinematic configuration when the fermion-antifermion pair of the anomaly 
diagram is collinear. However both in the case of the AVV diagram and in the conformal case, as 
evident from the expression of F\, it reappears as an extra contribution and is responsible for the trace 
anomaly. It is rather easy to show the pole dominance of the anomaly away from the conformal point 
(massive case) at high energy, since the non anomalous terms present in F% and F2 are subleading at 
large s. We are entitled to separate the pole contribution, which describes the non-local contribution 
to the trace anomaly, from the rest, and rewrite the F% form factor and effective action, respectively, 
as 

Fi = F lpole + F 1 (92) 

and 

S = S po i e + S. (93) 

The reminder (S) includes all the remaining contributions coming from the several form factors of the 
expansion, while the pole part gives 

S po ie = - / dAxdA y ( nh ^ ~ d^d u h^(x)) U- l y F aP {x)F^{y). (94) 

As we have just mentioned, it is not difficult to show that the anomaly pole in F\, in the general 
kinematical case (e.g. for off-shell photons and a massive fermion in the loop) decouples in the 
infrared (i.e. its residue vanishes) while it remains coupled in the massless on-shell limit. In other 
configurations (for any of the two photons off-shell) is also decoupled. This behaviour is in perfect 
analogy with the chiral case [20J. 

5.2 Massive and massless contributions to anomalous Ward identities and the 
trace anomaly 

Anomalous effects are associated with massless fermions, and for this reason, when we analyze the 
contribution to the anomaly for a massive correlator, we need to justify the distinction between 
massless and massive contributions. The latter contribute to the anomalous Ward identity, in our 
approach, via terms of 0(m 2 / s 2 ), where s = k 2 is the virtuality of the graviton vertex. At nonzero 
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momentum transfer (k 7^ 0) the second term on the right-hand side of Eq. ()85[) is interpreted as an 
anomalous contribution, proportional to an asymptotic (3 function (/3 as ) of the theory, coming from 
the residue of the anomaly pole which appears in the form factor F±. While the appearance of the 
asymptotic f3 function of the theory (which coincides with the (3 function of the MS scheme) is expected 
at large s, where all the remaining scales of the theory (si,S2,m) can be dropped, corrections to the 
asymptotic description in the ultraviolet (UV) are expected. At the same time, in the far infrared 
(IR) region, below the fermion mass, the anomalous contribution should approach zero in a certain 
fashion, which will be specified below. 

A complete quantitative understanding of this point for a general kinematics (e.g. for s / 0) 
remains, in a way, an open issue, but much more can be said for the simpler case of zero momentum 
transfer, where a consistent pattern of separation between massless and massive contributions to the 
correlator emerge in the UV region. In this case the virtuality of the two photons and the fermion mass 
m (plus a renormalization scale /U or M) are the scales which appear in the renormalized perturbative 
expansion. Related analysis have been presented in [27j and [TB] and our conclusions do not differ 
from these previous investigations. We summarize the main points. 

Respect to the case of the chiral anomaly, the trace anomaly is connected with the regularization 
procedure involved in the computation of the diagrams. In our analysis we have used dimensional 
regularization (DR) and we have imposed conservation of the vector currents, the symmetry require- 
ments on the correlator and the conservation of the energy momentum tensor. As we move from 4 to 
d spacetime dimensions (before that we renormalize the theory), the anomaly pole term appears quite 
naturally in the expression of the correlator. This is not surprising, since QED in d 7^ 4 dimensions is 
not even classically conformal invariant and the trace of the energy momentum tensor in the classical 
theory involves both a F 2 term (~ (d — 4)F 2 ) beside, for a massive correlator, a tpip contribution. 
Let's summarize the basic features concerning the renormalization property of the correlator as they 
emerge from our direct computation. 

1) The anomalous Ward identity obtained by tracing the correlator (r^ a ^) with involves 
only the Fi and F2 form factors in the massive case; in the massless case the scale breaking appears 
uniquely due to F\ via the term e 3 /(127r 2 )it Q ^(p, q), as pointed out before. The finiteness of the two 
form factors involved in the trace of the correlator is indeed evident. 2) The residue of the pole term 
(e 3 /(12-7r 2 )) in F\ is affected by the renormalization of the entire correlator (the form factor F\% is the 
only one requiring renormalization) only by the re-definition of the bare coupling (e 2 ) in terms of the 
renormalized coupling (e 2 ^) through the renormalization factor Z3. At this point, the interpretation 
of the residue at the pole as a contribution proportional to the f3 function of the theory is, in a way, 
ambiguous [28], since the j3 function is related to a given renormalization scheme. We stress once more 
that Eq. (|85p does not involve a renormalization scheme - which at this point has not yet been defined 
- but just a regularization. We have regulated the infinities of the theory but we have not specified a 
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subtraction of the infinities. For this reason, the substitution 



(e 3 /(12^ 2 )) -> 2/3 as (e)/e (95) 

which attributes the mass- independent term in F\ to a specific f3 function, the asymptotic one (/3 as ), 
as we are going to elaborate below, requires some clarification. 

To fully appreciate this point, it is convenient to go back to the unrenormalized Ward identity (|42p 
and differentiate it with respect to the momentum q and then set p = —q (k = 0) by going to zero 
momentum transfer. One obtains the derivative Ward identity 

9^™Hp, -p) = 2p 2 ^(p 2 )(p 2 g^ - pV). (96) 

The appearance of the derivative of the scalar self-energy of the photon on the right-had side of the 
previous equation is particularly illuminating, since it allows to relate this expression to a particular 
/3 function of the theory, which is not the asymptotic f3 as considered in Eq. (195j) . This j3 function is 
useful for describing the IR running of the coupling. 

To illustrate this point we start from the expression of the scalar amplitude appearing in the 
photon self-energy in DR 

2 \ e 2 f I 7 f 1 , M Nl _ m 2 — p 2 x(l 



n <" • = 5? \ja - i - J, ^ - *> "* — J < 97 > 

whose renormalization at zero momentum gives 

U R (p 2 ,m) =U{p 2 ,m) -II(0,m) = f x(l - x) log m ~ P x ^~ x \ ( 98 ) 

Using this expression, we can easily compute 

9 dU odU B e 2 e 2 m 2 f 1 , x(l — x) . . 

2p 2 - n = 2p 2 -f = / dx 2 \ 1 (99) 

dp A dp A D7T Z 7T Z Jo m z — p z x\l — X) 

Notice that this result does not depend on the renormalization scheme due to the presence of the 
derivative respect to p 2 . Notice also that the j3 function of the theory evaluated in the zero momentum 
subtraction scheme is exactly given by the right-hand side of the previous expression 



(where f3(e 2 ,p 2 ) = 2e[3(e,p 2 )), but this result does not hold, generically, in any scheme. The 
identification of anomalous (massless) effects in the theory, as exemplified by these simpler Ward 
identity, should then be obtained by extracting the appropriate (3 function of the theory, whose running 
should be driven by the effective massless degrees of freedoms (fermions, in our case) at the relevant 
observation scale (p 2 ). 
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Clearly, in the case of Eq. (jlOOp all the mass contributions have been absorbed into the very 
definition of the (3 function. Notice that if p 2 <C m 2 this (3 function, after a rearrangement gives 

^)-4/' fe pod 



e2 71-2 Jo m2 ~ P 2 x(l — x) 

and therefore it vanishes as j3 ~ 0(p 2 /m 2 ) for p 2 — >• 0. Equivalently, by taking the m — > oo limit 
we recover the expected decoupling of the fermion (due to a vanishing f3) since we are probing the 
correlator at a scale (p 2 ) which is not sufficient to resolve the contribution of the fermion loop. On 
the contrary, as p 2 — > oo, with m fixed, the running of the f3 function is the usual asymptotic one 
fiasie 2 ) ~ e 4 /(67r 2 ) modified by corrections 0(m 2 /p 2 ). The UV limit is characterized by the same 
running typical of the massless case, as expected. 

Notice that the right-hand side of Eq. (|96p . as we have already remarked, does not depend on the 
renormalization scheme, while the f3 function does and Eq. (jlOOp should be understood as a definition. 
For this reason, {3(e 2 ,p 2 ) correctly describes the IR running of the coupling as p 2 <C m 2 , and in this 
case it is obvious that massless anomalous effects of scale breaking are not present in this specific 
limit. 

In the case of regularization scheme different from zero momentum subtraction, there are some 
differences which should be taken into consideration. For instance, in a mass-dependent scheme one 
subtracts the value of the graph at a Euclidean momentum point p 2 = —M 2 , redefining the scalar 
self-energy as 



U R (p 2 ,m,M) = U(p 2 ,m)-U(p 2 = -M 2 ,m) = ^ ^ ( " ^ m * ^ ' 

2tt 2 



/ dx x(l — x) log 
Jo 



m 2 + M 2 x(l 

which gives, respect to the previous (M = 0) scheme, a /3 function now of the form 

e d e 2 f 1 m 2 —p 2 x(l — x) 



(102) 



/3(e) = dxx{l-x)\o g 



m 2 + M 2 rc(l 



e 3 f 1 , . , M 2 x(l-x) , . 

—Z dxx{l-x) V 2 ; 103 

2tt z Jo m z +M z x(l — x) 



For large values of M, this (3 function describes the usual UV running since 



e 3 f 1 , , , „, , e 3 



^ ~ 2^ / dXx{l ~ X) = /3(6)as = 12^2 • (104) 

J 

In this second scheme, the (regularization independent) right-hand side of Eq. (|96p can be interpreted 
as due to an anomalous contribution coming from the pole plus some explicit mass corrections, as 
obvious from the first and second term of (|99p . We conclude with some considerations on a third 
(mass-independent) scheme. 

In the MS scheme, the renormalization of the photon self-energy is performed via the subtraction 



■2 



U R (p 2 ,m, n) = U(p z , m, n) - ( - + 7 - log 4tt ) (105) 
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which gives directly an asymptotic /3 function since 




(106) 



It is clear, from these considerations, that a judicious definition of the /3 function allows a correct 
interpretation of the right-hand side of (|96p and (|99p . In the MS scheme, the breaking of scale 
invariance can be attributed to a UV running of the coupling (for p 2 3> m 2 ) plus mass corrections 
which are suppressed as 0(m 2 /p 2 ). Notice that in this case the renormalization scale (fi 2 ) should be 



sense, the extrapolation of the MS result to p 2 ~ fi 2 <C m 2 should be forbidden by the same criterion, 
since large logs of the relevant scales (log(m///)) would otherwise be generated. In the far infrared 
region p 2 <C m 2 the use of the same /3 function is indeed not appropriate, since the same scheme does 
not correctly describe the decoupling of the anomaly, which instead should occur, since there is no 
massless fermion in the theory. 

To conclude this discussion we just mention that the MS scheme can be used, obviously, both to 
describe the far IR and the far UV regions of the theory, with the condition that we are bound to 
choose a vanishing /3 function at p 2 <C m 2 and an asymptotic one for p 2 S> m 2 and assuring continuity 
of the gauge coupling across the fermion mass scale though the f3- function is discontinuous. This is the 
standard procedure followed in the MS scheme as, for instance, in QCD factorization, improved with 
the inclusion of threshold effects at the crossing scales (see for instance [29J, |30j) where the number 
of massless flavours change. 

5.3 The off-shell massless (TJJ) correlator 

Clearly, as we perform the massless limit on the amplitude, the residue of the same anomaly pole - 
identified above in the contribution F\ po i e - is still present, but will now be decoupled in the infrared. 

In the massless case the scalar functions Fi depend only on the kinematic invariants s,s\,S2 but 
we still retain the last entry of these functions and set it equal to for clarity, using the notation 
Fi = Fi(s; s±, S2, 0). These new functions are computed starting from the massive ones and letting 
m — > and Ao(m 2 ) — > 0, i.e. eliminating all the massless tadpoles generated in the zero fermion mass 
limit. 

The off-shell massless invariant amplitudes i^(s; si, S2, 0) are here given in terms of a new set of 
master integrals listed in Appendix [XJ We give here only the simplest invariant amplitudes, leaving 
the remaining ones to the appendix [Ej The anomaly pole is clearly present in F%, which is given by 



0(p 2 ), since we should not allow large logarithms to be present in the perturbative expansion. In this 



Fi(s; si, s 2 , 0) 




(107) 
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while 



F 2 (s; si, s 2 , 0) = 0. (108) 

The complete (TJJ) correlator is very complicated in this case as the long expressions of the form 
factors show, but a deeper analysis of its poles by computing the residue in s = can be useful to 
draw some conclusions. The single pieces of T fJ,ua/3 (s; si, s 2 , 0) indeed contribute as 

(109) 
(110) 
(111) 
(112) 



lim sF\(s; 


si, 


s 2 ,0)t^ = 


6 + fj.ua fS 




18 7T 2 h 


s=0' 


lim sF^ts; 


si, 


S2 ,0)tf^ = 


72 7T 2 3 ls 


=0' 


lim sF^(s; 

s->0 


si, 


a 2 ,0)tf^ = 


72 7T 2 5 ls 


=0' 


lim sFj(s; 

s^0 


si, 




P 2 

e .fiuap 1 

36vr 2 7 ls - 


=0' 



while F 2 is absent in the massless case. The residues of the si, s 2 , 0) not included in the 

equation above are all vanishing. Combining the results given above one can easily check that the 
entire correlator is completely free from anomaly poles as 

lim sT^ al3 (s; s lt a 2 , 0) = (113) 

in this rather general configuration. A similar result holds for the correlator responsible for the chiral 
anomaly and shows the decoupling of polar contributions in the infrared. 

5.4 The on-shell massive (TJJ) correlator 

A particular case of the (TJJ) correlator is represented by its on-shell version with a massive fermion in 
the loop. If we contract u a ^(p, q) and w a "{p, q) with the polarization tensors e a (p) and ep{q) requiring 
e a (p)p a = 0, £p(p)p^ = , the first tensor remains unchanged while w a ^(p,q) becomes w a/3 (p,q) = 
si S2g a P ■ This will be carefully taken into account when computing the s\ — > 0, s 2 — > limit of the 
product of the invariant amplitudes Fj with their corresponding tensors tj iual3 (i = 1, . . . , 13). 
The invariant amplitudes reported below describe Fj(s;0, 0, w?) whose tensors t^ va ^ are also finite 
and non-vanishing. They are 
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where the on-shell scalar integrals P(s, 0, 0, m 2 ) and Co(s, 0, 0, m 2 ) are computed in Appendix[Al here 
-^13 R denotes the renormalized amplitude, obtained by first removing the UV pole present in the 
photon self-energy by the usual renormalization of the photon wavefunction and then taking the on- 
shell limit. The remaining invariant amplitudes Fi(s, 0, 0, to 2 ) are zero or multiply vanishing tensors 
in this kinematical configuration so they do not contribute to the correlator. 

The limit from the massive on-shell form factors to the massless ones is clearer by looking at the series 
expansion of the scalar integrals around to = 

1 r. / s m 2 2 m 2 



2s 



C7oM,0,m 2 ) 
and from this we obtain for F[ 

F{(s,0,0,m 2 



log 



TO 



log 



+ C»(m d ) 



e 2 m 2 



1 



l r 
4 



log 



m- 



3vr 2 s 2 

where the notation F[ denotes the first form factor after the subtraction of the pole in 1/s. 
Using the results given above, the full massive on-shell amplitude is given by 
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so that the invariant amplitudes reduce from 13 to 3 and the three linear combinations of the tensors 
can be taken as a new basis 
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as previously done in the literature |22] . If we extract the residue of the full amplitude we realize that 
even though some functions Fi(s, 0, 0, to 2 ) have kinematical singularities in 1/s this polar structure is 
no longer present in the complete massive correlator 



limsT^ = 

s^0 
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showing that in the massive case the (TJJ) correlator exhibits no poles. In a following section we will 
comment on the interpretation of these massless poles exploiting the analogy with a similar situation 
encountered in the case of the gauge anomaly. 



6 The general effective action and its various limits 

In this section we present results for the correlator in various kinematical limits. We start from its 
expression in the on-shell massive case and then perform its expansion in 1/m which will be used 
in a next section to extract the corresponding effective action. As a final step we show the on-shell 
structure of the invariant amplitudes in the conformal limit. 

It is possible to identify from them the structure of the effective action in its most general form. If 
we denote by Si the contribution to the effective action due to each form factor Fi , then we can write 
it in the form 

Si = J d 4 xd 4 yd 4 zir"(z,x,y)h^(z)A a (x)Af,(y) J ^0^ e -^-~ z ^<y-^F t (k,p,q) (122) 

where k = p + q. We have introduced the operatorial version of the tensor structures tf ual3 , denoted 
by tj that will be characterized below. Defining 

* = « S! i' (123) 

and using the identity 

m z ,P x ,q y )S 4 (x - z)5\y - z) = J ^pV^*"* ^^F^p, q) (124) 

where formally Fi is the operatorial version of Fi, we can arrange the anomalous effective action also 
in the form 

S l = J d 4 xd 4 yd 4 zFi(k z ,p x ,q y ) [5 4 {x - z)5 4 (y - z)] (z,x,y)h^A a (x)A^y). (125) 
For instance we get 

i1™P(z,x,y)h^A a (x)Ap( y ) = l -{U z h{z)-d^d z u h^{z))F aP (x)F^(y) } (126) 

^ ,aP {z,x,y)h lw A a {x)Ap{y) = {U z h{z) - d;d*h^(z)) d^(x)d l/ F^(y), (127) 

i^(z,x,y)h^A a (x)Ap(y) = ^hT{z) {U x9lxv - 4d*c%) F a p(x)F a ?(y), (128) 

t^(z,x,y)h^A a (x)A p (y) = hT{z) (D x g^ - 40%%) d fl F^(x)d u F" x (y), (129) 

i^(z,x,y)h^A a (x)Ap(y) = ±h^(z) (D y g^ - 48»d») F a p{x)F^{y), (130) 

i^(z,x,y)h^A a (x)A p (y) = h^{z){U y g^-Adf l d«)d^{x)d v F v \y), (131) 
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t^(z,x,y)h^A a (x)Ap(y) = l -h^{z) (d^d^ - 2(3^ + 8*8%)) F a p(x)F^(y), (132) 
t^(z,x,y)h^A a (x)A p (y) = h^(z) j> A d^„ - 2(d»d^ + 3^)) d^(x)d„F»\y) (133) 

and similar expressions for the remaining tensor structures. However, the most useful forms of the 
effective action involve an expansion in the fermions mass, as in the l/m formulation (the Euler- 
Heisenberg form) or for small m. In this second case the non-local contributions obtained from the 
anomaly poles appear separated from the massive terms, showing the full-fledged implications of the 
anomaly. This second formulation allows a smooth massless limit, where the breaking of the conformal 
anomaly is entirely due to the massless fermion loops. 

In the l/m case, for on-shell gauge bosons, the result turns out to be particularly simple. We 
obtain 

7p 2 1 p 2 s 1 / 1 \ 

2160it z m z 3024 it m \m b J 

F 3 (s,0,0,m 2 ) = F 5 ( S ,0,0,m 2 ) = —^^+ f S — + O ( , (135) 

4320 ^to^ 60480 tt z m 4 \m D / 

F 7 (s,0,0,m 2 ) = -4F 3 (s,0,0,m 2 ) (136) 

F 13 , RS ,0,0,m 2 = ___ + ___ + _ , 137 
1440 tt z m z 20160 tt z m 4 \m D / 

which can be rearranged in terms of three independent tensor structures. Going to configuration 
space, the linearized expression of the contribution to the gravitational effective action due to the 
TJJ vertex, in this case, can be easily obtained in the form 



Stjj = J d 4 xd A yd 4 zT^(x,y, z)A a {x) A p (y) h^(z) 

J d 4 x {ah - Wh^) f 2 



4320 vr 2 m 2 

„2 



4320 vr 2 m 2 
lie 2 



d A x {uhF 2 - 8d^F al3 d v F a(3 h fMU + 4(d»d u F al3 )F al3 h^ 



+ 1440^/ <^DV. (138) 
which shows three independent contributions linear in the (weak) gravitational field. 



7 The massless (on-shell) (TJJ) correlator 

The non-local structure of the effective action, as we have pointed out in the previous sections, is not 
apparent within an expansion in 1 /m, nor this expansion has a smooth match with the massless case. 

The computation of the correlator F^ ual3 (s; 0, 0, 0) hides some subtleties in the massless fermion 
limit (with on-shell external photons), as the form factors Fi and the tensorial structures t{ both 
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contain the kinematical invariants si,S2- For this reason the limit of both factors (form factor and 
corresponding tensor structure) F{ t^ ua ^ has to be taken carefully, starting from the expression of the 
massless s%, S2, 0) listed in Appendix [E] and from the tensors t^ u contracted with the physical 
polarization tensors. In this case only few form factors survive and in particular 

Fi (a, 0,0,0) = 



F 3 (s, 0,0,0) 
F 7 (s, 0,0,0) 

F 13 p(s, 0,0,0) 



18ttV 
F 5 (s, 0,0,0) = - 
-4F 3 (s,0,0,0), 

-w f 121og 
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(139) 

(140) 
(141) 

(142) 



and hence the whole correlator with two onshell photons on the external lines is 
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where t t j LVa ^ are the tensors defined in eqs. (|118H120|) . 

The study of the singularities in 1/s for this correlator requires a different analysis for F± and 
the remaining form factors, as explicitly shown in eq. 11431 where F\ has been kept aside from the 
others, even if it is proportional to F%. Indeed iq is the only form factor multiplying a non zero trace 



tensor, t 



fj.va/3 
1 J 



and responsible for the trace anomaly. If we take the residue of the onshell correlator 
for physical polarizations of the photons in the final state we see how the 4 form factors and their 
tensors combine in such a way that the result is different from zero as 

„2 



lim sT^{s; 0,0,0) 



12tt 2 



■f q Q (p fl p u + q" q v 



(144) 



where clearly each singular part in 1/s present in iq , F 3 , F§ , F7 added up and the logarithmic behaviour 
in s of -F13 has been regulated by the factor s in front when taking the limit. The result shows that 
the pole, in this case, is coupled in the IR, as shown by the dispersive analysis. 



8 Conclusions 

We have presented a computation of the TJJ correlator, responsible for the appearance of gauge 
contributions to the conformal anomaly in the effective action of gravity. We have used our results to 
present the general form of the gauge contributions to this action, in the limit of a weak gravitational 
field. One interesting feature of this correlator is the presence of an anomaly pole |16j . 

Usually anomaly poles are interpreted as affecting the infrared region of the correlator and appear 
only in one special kinematical configuration, which requires massless fermions in the loop and on-shell 
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conditions for the external gauge lines. In general, however, the anomaly pole affects the UV region 
even if it is not coupled in the infrared. This surprising feature of the anomaly is present both in the 
case of the chiral anomaly [20] and in the conformal anomaly. Here we have extracted explicitly this 
behaviour by a general analysis of the correlator, extending our previous study of the chiral gauge 
anomaly. 

Indeed anomaly poles are the most interesting feature, at perturbative level, of the anomaly, being it 
conformal or chiral, and are described by mixed diagrams involving either a scalar (gravitational case) 
[16J or a pseudoscalar (chiral case) (20J, [3Tj . The connection between the infrared and the ultraviolet, 
signalled by the presence of these contributions, should not be too surprising in an anomalous context. 
The pole-like behaviour of an anomalous correlator is usually "captured" by a variational solution of 
a given anomaly equation, which implicitly assumes the presence of a pole term in the integrated 
functional [32]. By rediscovering the pole in perturbation theory, obviously, one can clearly conclude 
that variational solutions of the anomaly equations are indeed correct, although they miss homogeneous 
solutions to the Ward identity, that indeed must necessarily be identified by an off-shell perturbative 
analysis of the correlators. This is the approach followed here and in [20J. 

We have also seen that the identification of the massless anomaly pole allows to provide a "mixed" 
formulation of the effective action in which the pole is isolated from the remaining mass terms, ex- 
tracted in the S po i e part of the anomalous action, which could be used for further studies. We have 
also emphasized that a typical 1/m expansion of the anomalous effective action fails to convey fully 
the presence of scaleless contributions. 

There are various applications of our analysis which can be of interest for further studies. The 
first concerns the possible implication of these types of effective actions in cosmology, especially in 
inflationary scenarios where the coupling of gravity to matter via gauge interactions and the conformal 
anomaly plays an interesting phenomenological role. As we have stressed in the introduction, the 
local description of an anomalous effective action involves additional degrees of freedom which can 
be identified in the case of the gauge anomaly [3T] as well as in the conformal case In [IB] 

the authors describe the role of the corresponding scalar degrees of freedom in the effective action 
emphasizing their meaning as possible composite. In the case of an anomalous gauge theory the 
derivative coupling of the anomaly pole to the anomalous gauge current indicates that its nature is 
that of a quasi Nambu-Goldstone mode. If the parallel with chiral theories holds, one should be able 
to characterize the physical property of this state, including its BE condensation under the action of 
gravity. Similar features could be shared by the scalar state(s) described by the conformal anomaly. 
Of particular interest are the extensions of these analysis to the case of supersymmetric theories, in 
particular to N = 1 superconformal theories, where the R-current, the supersymmetry current and the 
energy momentum tensor belong to the same supermultiplet, as are their corresponding anomalies. 
Clearly our computation is the first step in this direction, and can be extended with the inclusion 
of other types of fields in the perturbative expansion, reaching, as a starting point, all the relevant 
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fields of the Standard Model. In general, one could also use our approach to come with a complete 
description of the interplay between supersymmetry and the conformal anomaly, acting as a mediator 
of the gravitational interaction, which is of phenomenological interest. Finally, we mention that our 
analysis could be useful in order to test, in a specific realization perturbative results of conformal 
field theories in four spacetime dimensions. These studies of the TJJ correlator have been performed 
on rather general grounds, using conformal invariance as a unique assumption in order to infer the 
structure of the operator product expansion [33J. We hope to return to these points in the near future. 
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A Appendix. Definitions and conventions for the scalar integrals 

We collect in this appendix all the scalar integrals involved in this computation. To set all our 
conventions, we start with the definition of the one-point function, or massive tadpole Ao(m 2 ), the 
massive bubble Bo(s,m 2 ) and the massive three-point function Cq(s, 8%, S2, m 2 ) 
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where S3 = s and in the last equation i = 1, 2, 3 and j, k ^ i. 

The one-point and two-point functions written before in n = 4 — 2 e are divergent in dimensional 
regularization with the singular parts given by 
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with 
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We use two finite combinations of scalar functions given by 
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The scalar integrals Cq{s, 0, 0, m 2 ) and 2?(s,0, 0,m 2 ) are the {si — > 0, S2 — >■ 0} limits of the generic 
functions Cq(s, s±, S2,m 2 ) and £>i(s, s\,m 2 ) 
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The master integrals denoted by Bq(s,Q), T>i(s,Si,0) (i = 1,2) and Cq(s, s±, s%, 0) are consistently 
redefined for m = (and s < 0) as 
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where /U is the renormalization scale and the function <f>(x,y) is defined as 

if y 1 + 7r 2 
= Ti 2[Li 2 (-px) +Li 2 (-/9y)j +ln-ln— h ln(px) ln(py) + — 



with 



A(x,y) = VA, 
p(x,y) = 2(1 - x-y + A)- 1 , 



A = (1 — x — y) 2 — Axy, 

si s 2 
x = — , y = — . 



(155) 
(156) 
(157) 

(158) 



(159) 
(160) 



s s 

The singularities in 1/e and the dependence on the renormalization scale /i cancel out when taking into 
account the difference of two functions Bq, so that the TVs are well-defined; the three-point master 
integral is convergent. 



30 



B Appendix. Alternative conditions on the correlator in the mass- 
less case 



As we have mentioned, one can follow an entirely different approach in order to fix the expression 
of the correlator. This is based on the requirement that the trace anomaly satisfies a well known 
operatorial relation which is imposed on the matrix elements at nonzero momentum. Specifically we 
proceed as follows, and illustrate this point in the massless limit. We impose the value of the trace 
anomaly as a defining condition on the whole amplitude, so that the (new) request c') will be 

c') the non-zero anomaly trace in the massless limit. 

As the first two conditions a) and b), respectively the {fi v} symmetry and the vector current 
conservation, remain the same as before, we continue illustrating the modifications due to this approach 
from this point on. The third condition is given by 

g^T^(p,q) = cu^(p,q), (161) 

where c is related to the usual QED f3- function as c = — A The resulting system is 

' -A 7 + 2A 9 -A 12 = 0, 

Eq. MB) < (162) 
c + 4 A 37 + 4 A A2 + A 4 p ■ p - 2 A 6 p ■ p + 2 An p ■ q + 2 A u q • q + A 16 q ■ q = 0, 

whose solutions for A^i and A37 read as 

A41 = ^-(A 7 -2A 9 + A 12 ) (163) 

^37 + ^42 = ~[-c-A A p-p + 2A 6 p-p-2A 11 p-q-2A u q-q-A 1(i q-q)}. (164) 

As seen from the last equation, the second solution returns the sum of two UV divergent amplitudes, 
^37 and A^2- However, an explicit computation shows that in the explicit mapping between the two 
sets of Ai and Fi these two amplitudes appear in such a way that their divergences cancel. Therefore, 
reinserting the expressions of An and A37 extracted from Eq. (I164D into the expression of T' iua ^(p, q) 
one finds another mapping between the form factors Ai and Fi, as previously done in Eqs. (I72ll84jl 

F ' " STP (165) 

F 2 = 0, (166) 

Fi = iF? <168) 
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^6 = j^, (170) 
Aq-q 



F c , A n (A 9 p-p + A u p-q)q-q p ■ p(A 6 p ■ q + A 9 q ■ H) 

7 6k-k + 2 + 2p-q> + 2p-q 2 ' [ ' 



F„ = -^ + A,,^\, (173) 
p • q p • q z 

F 10 = V^ + ^±, (174) 
p ■ q z p ■ q 

Fn = (175) 
2q ■ q 2p ■ p 

Fn = ^ + ^, (176) 
2p ■ q 2p ■ p 

F 13 = — — — \2Anp ■ q 2 + cp ■ q + 4 A i2 p ■ q + A 4 p ■ pp ■ q + 2A e p ■ pp ■ q 
Ap ■ q 1 

+ 2 A u q ■ qp ■ q + A^q • qp ■ q + 4A 9 p • pq ■ q\; (177) 

This new mapping leaves the invariant amplitudes from Fg to F\2 the same as before, so the condition 
c), i.e. the WI derived from Eq. (|36p and c') are perfectly equivalent in determining these 4 form 
factors. 



C Appendix. Comparison with the parametric approach and nu- 
merical checks 

The parametric approach of [16] allows, by combining the denominators of the various tensor am- 
plitudes, to give parametric expressions for the form factors F^ starting from a set scalar parametric 
integrals. Our results correspond to an explicit computation of these integrals. We will not give each 
integral separately, since they are rather lengthy. The mapping between the F^s in the parametric 
form and our expressions allow to perform numerical checks of our result. We have perfect agreement 
between the parametric forms derived in [16j . computed numerically, and our explicit expressions in 
all the euclidean regions of the external momenta. We briefly clarify this point. 

Explicit formulae for all twelve finite coefficient functions may be given in the Feynman parame- 
terized form, 

C 1 (k 2 ;p 2 ,q 2 ) = ^- 2 fdx r X dy— — 2 C -£^L— — (178) 

4vr z Jo Jo P x (l ~~ x ) + Q 2/(1 ~~ V) + 2 X D P " 1 + m 

where the polynomials (h(x, y) for i = 1, ... ,12 are listed in Table [3l 



Fi = C 7 + Cs + C 9 + p*_ + C3 + C8 _ Cg) + ^_ { _ Cj + Cs + Cw _ Ci2) j (1?9) 
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3 


Cj = coefficient of 


Cj(x,y) 


1 




-4x 2 (l -x)(l -2x) 
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(p^q u + q^p u )p a pP 


—x(l — x)(l — 4x + 8xy) + xy 


3 
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p^p u p a qP 
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q fJ/ q u p a qP 
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{p»q u + q^p u )q a pP 


-2xy(l-2x)(l-2y) 
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q^q v q a p fi 


2xy{\ - 2yf 


10 


p^p v q a q^ 


2y(l - 2x)(l - x - y + 2xy) 


11 


{p»q u + qVp»)q a qP 


-2/(1 - 2/)(l - 4y + 8xy) + xy 


12 


q^q u q a ql 3 


-Ay 2 (1 - 2y)(l - y) 



Table 3: The twelve tensors with four free indices (ixvafi) on p, q used in ref. |16| ) for the construction 



of the form factors At each coefficient functions Cj(k ,p ,q ) correspond a polynomial Cj in the 
Feynman parameterized form, as given in Eq. (|178p . 
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C*3 _ Cl2 
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Finally, numerical checks on F13 are performed on the UV convergent amplitude 

U R (p 2 )+U R (g 2 ) p 2 g 2 p 4 C 4 + g 4 C 6 p-g , 
-T13 = 1 ^5 H : 1 r - (,^2 + 03 + Oio + /On, 



4p • g 



+^ (2C 2 + C 4 + 2C 5 + C 10 ) + j (C 3 + 2C 5 + C 6 + 2Cn) 



(191) 



where the scalar two-point functions have been renormalized by the subtraction of the UV 1/e pole. 

D Appendix. The massive invariant amplitudes 

The off-shell massive form factors Fi, with 

• s / si^O s 2 ^0 m^O 
and with 7 = s — si — s 2 , a = s 2 — 2(si + s 2 ) s + (s\ — s 2 ) 2 are given by E] 
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1447T 2 scr 3 

„3\ „3 



s 6 - 3 (si - 4s 2 ) s 5 + 6 (3si - 7s 2 ) s 2 s 4 
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(194) 



F 4 (s; si, s 2 , 



F 6 (s; s 2 , si, m 



e 2 m 2 



3s 3 - 2 (2si + 3s 2 ) s 2 + (-sj + 6s 2 si + 3s 2 .) s + 2s x (s x - s 2 ) 2 

3s 3 -2(2si+3s 2 )s 2 + 



67r 2 sa 2 si 
e 2 

+ 10 2 2 [B (s, to 2 ) to 2 - *4 (m 2 )] 

127T z S(J z Sl 



35 



(-s? + 4s 2 s x + 3s 2 .) s + 2si (si - s 2 ) 2 
4 + (6s + llsi) s\ - (14s 2 + sis + 5s?) s 3 , + (16s 3 - 35sis 2 + 46s?s - 15s?) 



247r 2 cr 3 si 

- (s - si) 2 (9s 2 - llsis - 6s 2 ) s 2 + 2 (s - s x ) 4 (s + 2si) 



— e 2 2? 2 (s, s 2 , to 



m 



^ - 2 (2s + si) (s - si) 4 + (3s 2 - 43sis + 2s?) s 2 (s - si) : 



127T 2 S<7 3 Sl 

+ (9s + 4si) s\ + (-23s 2 + 29sis - 10s?) s\ + (l5s 3 + 2sis 2 + 5s?s + 6s?) sfj 

~ A8 J a 4 Si ( 3 ( s + s i) (^ - si) 6 - 4 (4s 2 - 14sis - 5s 2 ) s 2 (s - si) 4 
+ (35s 3 - 119sis 2 + 169s?s - 13s?) s 2 . (s - si) 2 + (s - 3si) s G 2 - 8 (s 2 + 9sis + 7s?) s\ 
+ (25s 3 + 159sis 2 - 197s?s + 157s?) sj + 4 (-10s 4 + 21s?s 2 + 28s?s - 27s 4 ) s^ 



+ e 2, Di(s, Si, m 



m 



2s 5 + (15si - 8s 2 ) s 4 + (-53s? - 5s 2 s x + 12s 2 .) s 3 



127T 2 S(T 3 Si 

+ (49s? + 46s 2 s? - 33s^si - 8s^) s 2 - (s x - s 2 ) (9s? + 52s 2 s? + 23s^si + 2s 3 ,) s 
-2si (si - s 2 ) 3 (2si + s 2 ) ^ + 48 J V4 f + 4 (6s + llsi) s 5 2 + (-87s 2 + 106sis - 91s : 
+4 (22s 3 - 69sis 2 + 40s?s + s?) sf, + (s - s x ) (3s 3 - 29s x s 2 + 209s?s - 79s?) s 2 2 
-8 (s - si) 3 (6s 2 - 13sis - 4s?) s 2 + (s - si) 5 (19s + 5si) 
4 r a(3s + 2si) 



+ e C (s,si,s 2 ,m ) 



TO, 
67T 2 CJ 2 



+ 18s 2 



m 



24tt 2 ct 3 



a 



s Sl 

+2si (si + s 2 )) + 12 (3s 2 - 3 (22si + 7s 2 ) s + 8l (3s x - 17s 2 )) a 
+720 s si ((s - si) 2 - 2 (s + si) s 2 ) 
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-79s 2 2 s 2 - 563s|s - 36s^)s? 
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where as previously done the master integrals are collected in Appendix These expressions have 
been analyzed in the text in various kinematical limits to show the appearance of anomaly poles and 
of all the other poles in the off-shell formulation. 

Notice that contains two vacuum polarization diagrams with different momenta on the external 
lines and has been renormalized by a subtraction at zero momentum 
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E Appendix. The massless invariant amplitudes 

We present here the expressions of the invariant amplitudes in the massless limit. We obtain 
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e si 



■ Di(a, si, 0) (s - si) 6 + 2 (14s + llsi) s 2 (s - si) 4 - (23s 2 - 214sis + 19s 2 ) s 2 . (s - si) ; 



48vr 2 a 4 

- 21s^ + 2 (5si - 2s) s| + (107s 2 - 318sis + 71s 2 ) + 8 (-lis 3 + 18sis 2 + 17s?s - 8s?) s^ 



e s 2 



48vr 2 



(7 



P 2 (s, s 2 , 0) \s\ - 2 (s - 14si) s| + (s 2 + 120sis - 37s?) s\ 



-4 (s 3 + 49sis 2 - 69s?s + 13s?) s\ + (s - si) (lis 3 - 69sis 2 + 309s?s - 83s?) s\ 
-2 (s - si) 3 (5s 2 - 49sis - 4s?) s 2 + 3 (s - s x ) 5 (s + 5s x ) 



16vr 2 



Co(s,si,s 2 ,0) 



2sis 2 



[2s 6 + 3 (s 2 - 3si) s 5 + (15s? + 6s 2 si - 13s|) s 4 



+2 (-5s? - 19s 2 s? + 29s^si + s 3 ,) s 3 + 12s 2 (4s? - 4s 2 s? - 3s^si + s 3 ,) s 2 
+ (si - s 2 ) 2 (3s? - 15s 2 s? - 31s 2 si - 5s 2 ) s — (si — s 2 ) 4 (si + s 2 ) 2 ] 
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F 4 (s; si, s 2 , 0) 



F 6 (s; s 2 , si, 0) 



., :j - s^ + (6s + llsi) s^ - (14s 2 + sis + 5s?) s\ + (16s 3 - 35s x s 2 + 46s?s - 15s?) s\ 

- (s - si) 2 (9s 2 - llsis - 6s?) s 2 + 2 (s - si) 4 (s + 2s x ) 

e 2 T 1 / 

^P 2 (s, s 2 , 0) - (3 (s + si) (s - si) 6 - 4 (4s 2 - 14sis - 5s?) s 2 (s - si) 4 

+ (35s 3 - 119sis 2 + 169s?s - 13s?) s\ (s - si) 2 + (s - 3si) s\ - 8 (s 2 + 9sis + 7s?) s\ 

+ (25s 3 + 159sis 2 - 197s?s + 157s?) s\ + 4 (-10s 4 + 21s?s 2 + 28s?s - 27s 4 ) si 

Pi(s,si,0) 



16vr 2 



^-j ( s^ + 4 (6s + llsi) s^ + (-87s 2 + 106sis - 91s?) s\ 



16vr 2 



+4 (22s 3 - 69sis 2 + 40s?s + s?) s\ + (s - si) (3s 3 - 29s x s 2 + 209s?s - 79s?) s\ 
-8 (s - si) 3 (6s 2 - 13sis - 4s?) s 2 + (s - si) 5 (19s + 5si) 

C (s,si,s 2 ,0) 



^ ( - 2s (s - si) 6 - 2 (s 2 + 7sis + 2s?) s 2 (s - s x ) 



+2 (7s 3 + 6sis 2 + lls?s - 4s?) s\ + 12 (2s 3 - 3s x s 2 - 2s?s + s?) s 2 . (s - s x ) : 
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-4 (s + si) sSj + 6 (s 2 - 5s lS + 2s 2 ) s\ - 4 (9s 4 - 25sis 3 + 33s?s 2 - 15s?s + 2s?) s 3 . 



F 7 (s; si, s 2 , 0) 

6si 



72vr 2 



340s s? 



[2(s + s 1 )s 2 -(s-s 1 ) 2 ) 



+ — ^ (-13s 2 + 166sis - 13s? + 39 (s + si) s 2 ) + 



3(-s + 27si + s 2 ) 



2 9s 6 

+ - + -*-- 



167T 2 



Co(s,si,s 2 ,0) <^ - 



4sis 2 



7 2 a 4 



cj s 7^ 7 

2 i oo„_„ /i „2\ „6 i n„ fr,p,2 on„.„ , oi„2 N \„5 



16vr 2 



Pi(s,si,0) 



(-9s 2 + 22s 2 s - 4s 2 .) s? + 2s (20s 2 - 39s 2 s + 21a|) s? 

+ (s - s 2 ) 6 s 2 (4s + s 2 ) + 2s (s - s 2 ) 4 (2s 2 + 5s 2 s + lis 2 .) si 
+ (-65s 4 + 96s 2 s 3 + 33s^s 2 - 62s^s + 6sf) s? 
+ 2s (27s 4 - 22s 2 s 3 - 108s^s 2 + W2s 3 2 s - 31s|) s? 

- (s - s 2 ) 2 (23s 4 + 40s 2 s 3 - 105sis 2 - 34s^s + 4s 2 1 )s? + s? - 

2si 



2ss? 



3 7 2 a 4 



(317s 2 + 227sis + 64s 2 ) s% + (s x - s) 5 (-7s 2 + 39sis + 32s 2 ) s 2 

+ (397s 3 + 846sis 2 - 539s?s + 312s?) s%-{s- si) 3 (23s 3 + 114sis 2 + 463s?s - 16s?) s 2 . 

- (275s 4 + si (1181s 3 + si (3si (93s + 94si) - 1441s 2 ))) sj + (s- si) (103s 4 + 767sis 3 

- 79s?s 2 - 563s?s - 36s 4 )s^ - 23s^ + (133s + 4si) s 7 2 + (s x - s) 7 (2s + s x ) 



16vr 2 



P 2 (s,s 2 ,0) 



2s 2 



3 7 V 4 



- 23s? + (133s + 4s 2 ) s\ - (317s 2 + 227s 2 s + 64s 2 .) s? 
+ (s - s 2 ) 5 (7s 2 - 39s 2 s - 32s 2 .) s 1 -(s- s 2 ) 7 (2s + s 2 ) 

+ (397s 3 + 846s 2 s 2 - 539s|s + 312s 3 ,) s? - (s - s 2 ) 3 (23s 3 + 114s 2 s 2 + 463s|s - 16s 3 ;) s? 
- (275s 4 + s 2 (1181s 3 + s 2 (3s 2 (93s + 94s 2 ) - 1441s 2 ))) sf + (s - s 2 ) (103s 4 + 767s 2 s 3 

-79s^s 2 - 563s^s - 36sf)s? 
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F 8 (s; si, 82, 0) 



4s| + (14si - lis) si + 2 (s - si) (2s + 9s x ) s 3 . 



16vr 2 



D 2 (s,s 2 ,0) 



127r 2 7<T 3 

+2 (7s 3 - 43sis 2 + 33s?s - 9s?) s\ - 2 (s - si) 2 (8s 2 - 21s x s - 7s?) s 2 
+ (s-si) 4 (5s + 4 Sl ) 

- — j - hs 7 2 + 3 (s - 13si) s^ + (57s 2 - 128sis + 43s?) s% 
37a 4 



- 5s^ + 3 (s - 13si) s% + (57s 2 - 128sis + 43s?) s\ 
+ (-155s 3 + 567sis 2 - 341s?s + 121s?) sj 

+3 (55s 4 - 176sis 3 + 86s?s 2 + 56s?s - 53s?) s^ + 3 (s - s x ) 6 (s + si) 
(s - si) 2 (75s 3 + 103sis 2 - 311s? s - lis?) s^ + (s - s x ) 4 (7s 2 + 124sis + 25s?) s 2 j 
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16vr 2 



Di(a,ai,0) 



37(7 4 



34 + 5 (5si - 3s) si + (27s 2 + 24s x s + lis?) s\ - (l5s 3 + 339sis 2 

-289s 2 s + 159s?) si + (-15s 4 + 616sis 3 - 714s 2 s 2 + 168s?s + 121s 4 ) si 
+ (s - si) (27s 4 - 402sis 3 + 40s?s 2 + 298s?s - 43s 4 ) s 2 2 

- (s - si) 3 (15s 3 - 51s x s 2 - 245s?s - 39s?) s 2 + (s - si) 5 (3s 2 + 22s lS + 5s 2 ) | 



16vr 2 



C (s 



, si, s 2 , 0) I -^j (si + s 2 ) s 7 - 6 (s? - s 2 si + s 2 ) s 6 



+ 3 (si + s 2 ) (5s 2 - 12s 2 si + hs 2 2 ) s 5 + 2 (-10s 4 + 3s 2 s? + 54s^s 2 + 3s^si - 10s^) s 4 
+ (si + s 2 ) (15s 4 + 16s 2 s? - 126s^s 2 + 16sijsi + 15sf) s 3 
- 6 (s? + 5s 2 s? - s^s 4 - 18s^s? - s|s 2 + 5s^si + sf) s 2 
+ (si - S2) 2 («i + S2) (4 + 6s 2 s? + 34s|«? + 6s^si + sj) s 
+ 2si (si - s 2 ) 4 s 2 (si + s 2 ) 2 



F 9 (s; si, s 2 , 0) 



Fi (s; s 2 , si, 0) 



12vr 2 7 2 a 2 _ 



(s - si) 4 - 4 (4s + si) s 2 (s - si) ; 



16vr 2 



C (s,si,s 2 ,0) 



- 3s^ + 4 (si - 2s) si + 2 (13s 2 - 2s x s + s?) si 
^| (si + (2si - 3s) s 4 + 2 (s - si) (s + 2si) s 3 



+ 2 (s 3 - 7sis 2 + 3s?s - s?) s 2 - (s - 3si) (s - si) 2 (3s + si) s 2 + s (s - si) 



16vr 2 



P 2 (s,s 2 ,0) 



4s 2 



37 2 o"' 



2^3 



(4 ~ (35s + llsi) si + 30 (3s 2 + s?) si + 2 (-35s 3 + 17sis 

+ llsf s - 17s?) sf. + (s - si) 2 (5s 2 + 26sis + 17s 2 ) s 2 + 3 (s - si) 4 (3s - si 
4 



16vr 2 



Pi(s,si,0) 



3 7 2 



^3 f - s^ + (18s + llsi) si - 3 (21s 2 - 3sis + 10s 2 ) s| 
+ 2 (46s 3 - 37sis 2 + 2s?s + 17s?) s 3 , - (63s 4 - 82sis 3 + 2s?s + 17s 4 ) a| 
+3 (s - si) 3 (6s 2 + 7sis - sf) s 2 - s (s - si) 5 



Fn(s; si, s 2 , 0) = F i2 (s; s 2 , si, 0) 



12tt 2 (t 2 s 2 



2s 3 - (5si + 2s 2 ) s 2 + (4sf + 4s 2 si - 2s 2 2 ) s - (si - 2s 2 ) (s x - s 2 ) : 



16vr 2 



C (s,si,s 2 ,0) 



4s 

'^3 



s 4 + (si - 4s 2 ) s 3 - (si - s 2 ) (5si + 6s 2 ) s 2 
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(208) 
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+ (si + s 2 ) (3s? + 3s 2 si - 4s 2 ) s + (si — s 2 ) 2 s 2 (3si + s 2 )^ 



16vr 2 



P 2 (s,s 2 ,0) 



^3 ^-17s 4 + (26si + 48s 2 ) s 3 - 42s 2 ( Sl + s 2 ) s 2 



-2 (si - s 2 ) (5s 2 + 17s 2 si + 4s^) s + (si - 3s 2 ) (s x - s 2 ) ; 



16vr 2 



Pi(s,si,0) 



3d 3 s 2 



3s 5 - (10si + 9s 2 )s 4 



+ 2 (6s 2 + 26s 2 si + 3s^) s 3 - 6 (sf + 4s 2 s 2 + Us 2 2 si - sfj s 2 

+ (si - s 2 ) (s? - 19s 2 s 2 - 43s 2 ,s i + 9sl) s + (s x - 3s 2 ) (s x - s 2 ) 3 s 2 



Fi3,r(s; si, s 2 , 0) 



i[n Ji (si,o) + n fl ( S2 ,o)] + ^ 



S7 3s 
— + — + 1 



+ 



16vr 2 



C (s,si,s 2 ,0) 



<j 7 

4s 2 sis 2 (s 2 - 2 (si + s 2 ) s + sf + s|) 
7a 2 



16vr 2 



Pi(s,si,0) 



370"^ 



> - si) (5s 3 + sfs - 4s?) s 2 + (10s 2 + 5sis + 7s 2 ) si 



(s - si) 3 (s 2 + 2sis - sf) + (-10s 3 + 3sis 2 - 7s?) s 2 2 + sl- (5s + 4s x ) sj 



16vr 2 



P 2 (s,s 2 ,0) 



370- 



■ ^ (4s 3 + sis 2 + 7s?) si. - (8s 2 + 5sis + 7s 2 ) s^ 



(s - si) 5 + (s - 4si) (s + si) s 2 (s - si) 2 - s% + (5s + 4si) s\ 



(209) 



(210) 



as already noticed above for the case of the massive form factors the last one, i.e. i*i3,.R(s; s±, s 2 , 0), has 
been affected by the renormalization procedure for which the one-loop transverse photon propagator 
with a virtual pair of massless fermions is given by 



n fl (s,o) 



12 7T 2 



(211) 



where the dependence on the renormalization scale \i remains explicit. 



F Appendix. The asymptotic behavior of the off-shell massless 
(TJJ) correlator 

We present here the asymptotic expression of the form factor in the high energy limit. The leading 
contributions to the expansion in each expression come from the pole singularities (conformal or 
anomalous) except for F13 which has a constant asymptotic term. 

e 2 

Fi(s,si,s 2 ,0) = -j^r s , 



44 



F 2 (s,s 1 ,s 2 ,0) = 
F3(s,s 1 ,s 2 ,0) = 

F 4 (s, si, s 2 ,0) = 



F 7 (s,si,s 2 ,0) = 

Fg(s,si,s 2 ,0) = 

F 9 (s,s 1 ,s 2 ,0) = 

F n (s,s 1 ,s 2 ,0) = 

Fi 3 (s,s 1 ,s 2 ,0) = 



0, 



144tt 2 s 48ttV 

" 2 



si + 6s 2 + «i log (y) + 3s 2 log 



487T 2 SlS 



31og(^) !- ! 



+ 



487T 2 SlS 2 



2vr 2 si + 16si + 6s 2 + 19si log f — 

V s 



+ log (^) (9*1 + 8s 2 + 6si log (^) ) 



e 2 e 2 
+ 



36vr 2 s 24vr 2 s 2 



3 (si + 8 2 ) + 28i log (y) + 2s 2 log (^) 



+ 



24vr 2 s 2 



(— ) + 31og(— ) 10 



+ 



1 



12vr 2 s 2 
e 2 

247T 2 S 2 S 



(7) 



(7) 



31og(y ) + ! 



+ 



247T 2 S 2 S 2 



6.s 



1 + 2^ 2 s 2 + 12s 2 + 17s 2 log (^j 



+ log (^j (8 Sl + 9s 2 + 6s 2 log 
-l[U R ( Sl ,0) + U R (s 2 ,0)] + ^ 



12vr 2 s 

„2 



log( i) +1 „ g ( ? ) + | 



'£2 
S 



si + s 2 + 2si log + 2s 2 log (" 

2 (s 2 + (3 + vr 2 ) s 2 si + si) + s 2 (13ai + 6s 2 ) log 



24vr 2 s 2 

+si log (^) (6si + 13s 2 + 6s 2 log (^)) 



+ 



(212) 



G Appendix. The asymptotic behavior of the on-shell massive (TJJ) 
correlator 

This appendix contains the asymptotic expansion of the relevant on-shell massive form factors, that 
is their dominant contributions as s — > 00 with s > after taking into account the suitable analytic 
continuation. They result 



F 1 (s,0,0,m 2 ) = + 



e 2 m 2 



18vr 2 s 12vr 2 s 2 



4 - log z ( — J - 2 m log f — ) + 7T 2 



- O [ 4 



(213) 



F 3 (s,0,0,m 2 ) = F 5 (s,0,0,m 2 



e 2 m 2 



144vr 2 s 24vrV 



log 2 



m 



m 



+ 2in) log ( — ) + - - 
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F 7 M,0,m 2 ) 
F 13 ,i?M,0,m 2 ) 



-6in - 14 



-4F 3 (s,0,0,m 2 ), 



+ 0[-s ), 



(214) 
(215) 



144 7T 2 



+(2 + 2i^)log 



12 log ( ^- ) f:?5 1 2 / 7T 



2 4 

em 

47T 2 S 2 



+ 



e 2 m 2 



log 2 



+ 10 - 7T 2 + 2ivr 



-log 2 



m 



+ (2 -2m) log ( — j -\- K - 
s l \ s 



+2m - 3 



+ o 



(216) 
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